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Abstract. - A new experimental procedure is discussed, which aims at separating thermal 
from quantum behavior independently of the energy barrier distribution in small particle 
systems. Magnetization relaxation data measured between 60 mK and 5 K on a sample of 
nanoparticles is presented. The comparison between experimental data and numerical calcula- 
tions shows a clear departure from a standard thermal dynamics scenario, a result which was 
not obvious without using the new procedure presented here. 



The prediction that the magnetic moment of a single domain particle should flip by (quasi- 
macroscopic) quantum tunneling through the anisotropy barrier jjj || has motivated many 
experiments. Some attempts at a study of a single particle have been made [||, but most of 
the experiments have been carried out on a macroscopic number of nanoparticles dispersed 
in a non magnetic matrix, mainly by measuring the magnetic relaxation after a field change 
(so-called "viscosity" measurements) || . In such systems the sizes and the anisotropy constants 
of the particles are distributed, and an accurate knowledge of these distributions is out of 
reach, particularly for the very small barriers which are of interest when measuring the 
slow dynamics at low temperature. In addition, there may exist numerous small energy 
barriers due to surface defects [|| 0), making hazardous a simple correlation between the 
energy barriers and the particle sizes. Consequently, the experimental evidence of quantum 
tunneling of the magnetization (QTM) based on viscosity measurements in such systems 
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remains controversial (jl pi. In this letter, we present relaxation measurements on small isolated 
magnetic particles at low temperatures. We describe in the first part the usual analysis of the 
thermal variation of the magnetic viscosity, emphasizing its limits. In the second part, we 
discuss a new experimental procedure, which is likely to give much more reliable information 
on the thermal or quantum nature of the observed phenomena. 

Viscosity measurements and their limitations. - The sample consists of small ferrimagnetic 
particles of 7 — Fe20s (maghemite) which are dispersed in a silica matrix, with a volume 
fraction 4.3 10~ 4 . A Transmission Electron Microscopy study shows that their sizes can be 
fitted to a lognormal distribution with peak value do = 6.3 nm and standard deviation a = 0.25 
[|] . It may thus be considered a good example of isolated single domain particle system. The 
measurements were taken with a home made combination of an r.f. SQUID magnetometer and 
a dilution refrigerator fll|| . The sample is coupled to the mixing chamber through a thermal 
impedance which allows a sample temperature range of 35 mK to 7 K . A 62 Oe magnetic 
field is applied on the sample at an initial temperature of 6 K, which proved to produce a well 
defined initial state for subsequent measurements at lower temperatures; as a matter of fact, 
at 6 K, all the processes that one can measure at lower temperatures are rapidly reaching 
their equilibrium. This will appear clearly from the simple calculations below, and has been 
checked by other choices than 6 K 0. After field cooling the sample down to the desired 
temperature and waiting for thermal equilibrium^), the field is cut off and the time change 
of the magnetic moment is recorded by following the variation of the SQUID signal( 2 ). This 
method avoids heating problems related to the sample movement at low temperatures, and 
allows a better sensitivity. We have studied the variation of the moment as a function of the 
time t after the field cut-off. On one decade of time, between 100 s and 1000 s, the plot 
is roughly linear in log(i). Fig.l shows the average logarithmic slope (viscosity S) of M(t) 
between 100 s and 1000 s for temperatures from 5 K down to 60 mK. 

For decreasing temperatures, the measured relaxation rate first decreases, then flattens out 
and surprisingly increases back below 150 mK. The same result is obtained with a ten times 
smaller preparation field (the viscosity being proportionally reduced). 
For one isolated uniaxial particle of anisotropy barrier U, the relaxation time can be written 

where To is a microscopic attempt time of the order of 10 -8 s to 10 -12 s. T*(T) is an 
effective temperature, which equals T in the case of thermally activated dynamics. The 
first-order predicted effect of a crossover towards quantum dynamics is that T*(T) should 
become greater than T or even temperature independent for temperatures below a certain 
crossover temperature T cr Q, leading to faster fluctuations than expected from thermal 
dynamics. Thinking of the sample relaxation at temperature To as a sum of independent 
processes (in our case magnetization reversal of isolated particles), one may write the total 
magnetic moment M(t, To) as 

M(t,T ) = £ mi (U) P(U) exp ( - ^phjTj) dU (2) 

( 1 ) By measuring the paramagnetic component of the sample signal [g| , we have checked that the 
sample temperature accurately follows the thermometer temperature in the whole accessible range, 
with a time constant less than 60 s. 

( 2 ) The time spent before the field cut-off has been checked to be of no significant influence on the 
relaxation rate, in contrast with what has been observed in more concentrated systems, where dipolar 
interactions may yield aging effects |LC]] . 
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Fig. 1. - Thermal variation of the measured relaxation rate. A 62 Oe field was cut off at t = s. 
Very low T zoom : x = rate around 100 s; • = rate around 1000 s . 



where P(U) is the energy barrier distribution, and m,i(U) is the average initial magnetic 
moment of the objects of anisotropy barrier U. From eq.^| and using the usual step function 
approximation ]l2| , the logarithmic relaxation rate approximates to 

S = - dM }^ 0) si k B T*(T Q ) P(U C ) mi (U c ) , where U c = k B T*(T ) In - . (3) 
a In t ro 

U c is the barrier energy of the objects having their main contribution to the dynamics after time 
t at temperature To. The distribution P(U C ) mi(U c ) has no reason to have a weaker dependence 
on To than T*(To) itself. The crucial point is that the relevant part of this distribution is out 
of reach in such a system : at low T and under zero-field, those entities which contribute to the 
dynamics correspond to very small energy barriers, whose physical origin remains uncertain 
[U q|. Therefore, extracting any T*(To) (which is the quantity of interest to characterize 
a departure from thermal dynamics) directly from a measurement of S(T) is generally not 
justified, for this implies arbitrary hypotheses over P{U) m,{U) @, |, §. In the next section, 
we present a new experimental procedure for reliably distinguishing between thermal and 
quantum dynamics, a procedure which is almost insensitive to the distribution P{U) rrii(U). 

Disentangling thermal from guantum dynamics. The point is that the temperature 
dependence of the relaxation rate is much weaker in the quantum regime than in the thermal 
regime. Hence the following idea : the thermal part of the relaxation towards equilibrium 
should rapidly be exhausted by a pre-relaxation at a higher temperature, which in contrast 
should be of little influence on quantum processes. Such a procedure is sketched in fig. 2a. 
After field cooling the sample from 6 if to a temperature x.Tq higher than To, we cut off the 
field, wait for to = 200 s and then cool down the sample to To. Still considering independent 
relaxation processes, and neglecting as a first approximation a possible x-dependence of the 
initial moments, we write the relaxation of the total magnetic moment as 

M(t,T , X ) = P(U) m t{ U) exp ( - exp ( - -*^) dU . (4) 
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As a first step, we now give a qualitative interpretation of our procedure, using the same kind 
of approximation as in eq.||. The logarithmic slope S(Tq, x) of M around time 1000 s after the 
field cut-off roughly corresponds to the product of the standard S(Tq) with a damping factor 
due to the pre-relaxation at x.Tq. By dividing S(To,x) by S(Tq), P(U c ) nn(U c ) is eliminated 
(see eq.^J), and we obtain a quantity which we may call "residual memory ratio" (RMR), in 
the sense that it represents the memory of the initial state that the system has kept through 
the whole procedure. 

Beyond this only qualitative argument, we have performed a complete calculation of eq.[|, 
and checked that the variation of this ratio RMR(x) = S(T , x)/ S(T ) indeed contains a 
very strong information upon T*(T), while being remarkably insensitive to the P{U) rrii(U) 
distribution. We have calculated the time variation of M and its average logarithmic slope 
between log(f) = 3.0 and log(f) = 3.2 (this choice has been checked to be of no influence and 
is related to the experimental procedure described below). In fig. 2b is plotted the calculated 
RMR(x) for various hypotheses on T*(T) and P(U) m l (U), with r = 10~ 10 s . 
In the thermal case (fig. 2b curves (i)), RMR{x) is the same for all working temperatures To; in 
a plateau hypothesis for quantum dynamics (fig. 2b curves (ii)), RMR{x) shows a corresponding 
plateau, followed by a sharp decrease at x = T cr /T . 




Fig. 2. - a) Typical magnetothermal procedure (here x = 1.2). The field is cut off at time t = s. 
b) Calculated RMR(x) for (i) thermally activated dynamics and (ii) quantum dynamics (in the plateau 
hypothesis with an arbitrary choice of T cr = 1 K and To = 0.5 K). The insert shows T*(T) in each 
case. Three choices of P(U) rrii(U) are presented. 



In both cases, the curves are very weakly dependent on the To value; for instance, in the 
thermal case, we have computed that the x-value at which RMR{x) has decreased by 90% 
ranges from 1.15 for r = 10 _8 s to 1.20 for r = 10 _12 s. 

A crucial result is that, as is clear from fig. 2b, the calculated RMR(x) is nearly insensitive to 
the extremely broad choice of P(U) m,i(U) presented here (ranging from U~ 5 to U +5 ), while 
it clearly reflects the thermal or quantum nature of the dynamics. Thus, within our present 
description, thermal dynamics can be characterized by a sharp decrease of RMR at low x. We 
have assumed temperature independent barriers U; in that respect, interactions between the 
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particles must be negligible. This was our motivation for choosing a sample of highly diluted 
particles. Also, we do not consider the case of a singular P(U) rrii(U); for instance, in the 
case of a delta- function peaking at U ^> U c (see eq||), the thermal RMR(x) would exhibit a 
slower x-decrease, but this would imply a very sharp temperature dependence of the viscosity, 
in disagreement with our present result. 

We have applied the procedure sketched in fig. 2a to measurements on the same sample as in 
fig.l; the results are shown in fig. 3a. 




Fig. 3. - a) Measured RMR(x) from 3 K down to 60 mK (symbols). The solid line is the calculation 
for thermally activated dynamics (see fig. 2b (i)). b) Calculated RMR(x) for the T*(T) shown in the 
insert, for different choices of To. 



For To = 3 K and To = 2 K, RMR(x) is very close to what is expected in the thermal 
regime. For lower temperatures, the observed RMR(x) can no more be explained in terms of 
thermal dynamics : its variation with x becomes slower and slower as temperature decreases, 
in an intermediate fashion between the two extreme cases plotted in fig. 2b. The behavior 
clearly departs from the thermal case of fig.2b(i); however, when compared with the ideal 
plateau case of fig.2b(ii), the smooth decrease of RMR(x) is suggestive of a distribution of 
crossover temperatures T cr in the system, or of a possible influence of an x-dependence of the 
initial moments. 

Keeping the same hypotheses as for eq.4, we have computed RMR(x) for a simplistic T*(T) 
shape as displayed in fig. 3b. We have chosen this T*(T) shape in order to obtain a qualitative 
agreement with our present experimental data. This T*(T) would still be of about 0.5 K at our 
lowest measurement temperature (0.06 K); this is far beyond any experimental uncertainty on 
the temperature of our sample. The T*(T) shape, together with the measured S(T) variation 
of fig.l, yields through eq.|| an estimate of the distribution P(U) rrii(U). It increases towards 
small energy barriers, as suggested but not proven in |^]. This calls for a better knowledge of 
the physical origin of low energy barriers in such samples. 

As a conclusion, this new procedure brings significant information on the nature of the dynam- 
ics, almost independently of the magnetic moment distribution P{U) rrii(U). The method has 
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allowed us to evidence a significant departure from the expected thermal dynamics scenario in 
an assembly of small magnetic particles. In our opinion, this method should help evidencing 
and studying QTM in many-particle samples, and could be interestingly extended to the 
characterization of non-thermal behavior in other systems, such as depinning of vortices in 
superconductors or of Bloch walls in ferromagnets. 

*** 
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